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We have proposed a method for correcting the Kohn-Sham eigen energies in the density functional 
theory (DFT) based on the Koopmans’ theorem using Wannier functions. The method provides a 
general approach applicable for molecules and solids for electronic structure calculations. It does 
not have any adjustable parameters and the computational cost is at the DFT level. For solids, the 
calculated eigen energies agree well with the experiments for not only the band gaps, but also the 
energies of other valence and conduction bands. For molecules, the calculated eigen energies agree 
well with the experimental ionization potentials and electron affinities, and show better trends than 
the traditional Delta-self-consistent-field (ASCF) results. 
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Density functional theory (DFT) l] is the main work¬ 
ing horse for material simulations, especially for ground- 
state properties such as atomic structures and binding 
energies. However, it is well known that the DFT, in 
particular the Kohn-Sham eigen energy [2], significantly 
underestimates band gaps. This is related to the lack of 
derivative discontinuity in the exchange-correlation (XC) 
potential when the total number of electrons crosses an 
integer point 3,0]. Over the years, various methods have 
been developed to overcome this deficiency. One popular 
approach is the hybrid functional , which mixes the 
exact exchange with local/semilocal XC functionals. Al¬ 
though widely successful, these methods depend on the 
mixing parameters and are computationally more expen¬ 
sive than local/semilocal functionals such as the local 
density approximation (LDA) or generalized gradient ap¬ 
proximation. An even higher-level method is the GW and 
the related random-phase approximation (RPA) 1,1, 
but the high computational cost of RPA makes it only 
applicable to small molecules [To}- There exists another 
approach to correct the DFT and its related Kohn-Sham 
Hamiltonian: the Koopmans’ theorem [nun. Under 
the Koopmans’ theorem, the total energy E(n) as a func¬ 
tion of the continuous number of electrons n should be 
a straight line between two integer points (11, [19|. How¬ 
ever, the LDA E(n) curve is convex [2(J, |2l|. Accord¬ 
ing to the Janak’s theorem, for a ./V-electron system, 

evBM = n— ►JV—, e CBM = ^r\ „-»jv+ with q being 

the eigen energies of the valence band maximum (VBM) 
and conduction band minimum (CBM) [20l422j. so the 
convexity of the LDA energy leads to an underestima¬ 
tion of the band gap ecbm — evbm- One explicit way to 
enforce the straight-line condition is to modify the LDA 
total energy to 

E = U LDA + ^ Msi) , Sl G [0,1] (1) 

i 

with E t {si) = {Ei(N ± 1 ) - E(N)\ x si [Ei(N ± si) - 
E(N)} (+ for adding electrons to unoccupied orbitals, 


and — for removing electrons from occupied orbitals). 
Ei{N ± si) is the self-consistent-field LDA energy after 
adding/removing si electrons on the 4>i orbital. Note 
Ei(si) = 0 for si = 0 or 1. Taking the variational mini¬ 
mum of E with respect to %pi (here s; = \{4>i\4>i)\ 2 is the 
projection of ipi on </>;), we obtain a modified Kohn-Sham 
equation 




( 2 ) 


with A i = dE Qg Sl ^ jgf=o, providing corrections to eigen en¬ 
ergies. 

Using eigen orbitals as (pi, the eigen energies in Eq.([2]) 
work successfully for atoms and small molecules 1 44171 ], 
because the straight-line condition makes the VBM and 
CBM eigen energies equal to the LDA total-energy differ¬ 
ences E(N)—E(N—1) and E(N+1)—E(N) (the so-called 
ASCF method [23|), which have been shown to agree well 
with the experimental atomic/molecular ionization ener¬ 
gies and electron affinities [24]. However, the above pro¬ 
cedure fails for extended systems such as solids, because 
adding/removing electrons on an extended eigen state 
only gives an infinitesimal charge-density change locally, 
which leads to Ei(si) = 0 and A; = 0 (no correction). In 
Ref. [iil ]. the authors proposed adding/removing a finite 
number of electrons in the primitive cell, which improves 
the LDA band gaps. However, an empirical parameter is 
needed to adjust the amount of electrons to be changed. 
So far, there is no systematic and parameter-free method 
to use the Koopmans’ theorem for correcting band gaps 
from molecules to solids, and the applicability of this ap¬ 
proach is significantly limited. 

In this work, we propose a new approach that ex¬ 
tends the above correction to solids. The basic ansatz is 
that the energy curve Ei(N ± si) for removing/adding si 
fractional electrons on any single-particle orbital (f>i (not 
necessarily eigen states) in the valence-band/conduction- 
band subspace should be a straight line. We call it an 
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extended Koopmans’ theorem, as it releases the condition 
that <fii must be an eigen state. This ansatz is plausible 
as in the Hartree-Fock total-energy expression, when the 
orbital relaxations are ignored, Ei(N ± si) is linear with 
si regardless whether (f>i is an eigen state or not, because 
the nonlinear self-interaction term is cancelled out be¬ 
tween the Hartree energy and the exchange integral. In 
this sense, the straight-line condition removes the self¬ 
interaction energy for (pi, which may improve the band 
gap [26j. Of course, the correctness of the ansatz can 
only be judged by its final results. To introduce nonzero 
corrections, we need localized (pi. The Wannier func¬ 
tions (wf) usm are most localized orbitals within the 
valence-band and conduction-band subspaces. Hence we 
will use WF as (pi in Eq. m and ©• We will show that 
the resulting eigen energies are in excellent agreement 
with experiments. 


The DFT calculations are performed by the PEtot 
code [29j within the LDA. Norm-conserving pseudopo¬ 
tentials are used. The plane-wave energy cutoff and k- 
point mesh guarantee the convergence of the eigen ener¬ 
gies within 0.01 eV. For the LDA band gap calculation, 
we use the primitive cell of the compound with the lat¬ 
tice parameters and atomic coordinates fully optimized. 
The spin-orbit coupling (SOC) is considered. 


The WFs are calculated using the wannier90 code [30[. 
We construct the WFs for valence bands and conduction 
bands separately, by projecting the eigen states onto the 


atomic orbitals [3JJ, |32|. The atomic orbitals are chosen 
based on their predominance in the valence and conduc¬ 
tion bands. For example, for GaAs, the valence-band 
WFs are constructed by projections onto one As s and 
three As p orbitals; the conduction-band WFs are con¬ 
structed by projections onto one Ga s and three Ga p or¬ 
bitals. The resulting WFs are shown in the Supplemental 
Material Fig. SI, which are all highly localized. We have 
tested the projected WFs and maximally-localized WFs, 
and the resulting band-gap corrections differ by less than 
0.01 eV. 


The next task is to calculate A/. We find that the 
screening effect from other electrons plays an important 
role. If we remove the WF orbital from the charge density 
and perform non-self-consistent calculations, the band 
gap corrections can be overestimated by several eV. To 
include the screening effect, we must calculate Ei{N±si ) 
self-consistently. We perform this calculation with spin- 
polarization using a supercell equal to 4 x 4 x 4 times 
primary cell, to exclude the interactions between the WF 
and its images. To add/remove fractional WF (pi in the 
spin-up channel, we optimize all the spin-up states or¬ 
thogonal to (f>i and optimize the spin-down states in the 


conventional way. The energy Ei{N ± si) is 

Ei(N ± si) = — 77V 2 + Vnl I “Pj ) + &i( 4 >i\ — T^V 2 + Vnl\<Pi) 

i 

p{ r ) = + on\(pi \ 2 

3 

where Vnl is the nonlocal-potential operator, Vi on (r) is 
the ionic potential, and Eh X c is the conventional Hartree 
and LDA XC energy. For valence-band WFs, ai = 1 — sj 
and the summation of j is from 1 to N/ 2 — 1 (N is an even 
number) for the spin-up channel, and 1 to N/2 for the 
spin-down channel; for conduction-band WFs, cq = Si, 
and the summation is from 1 to N/2 for both spin-up 
and spin-down channels. The ipj in the spin-up chan¬ 
nel is required to be orthogonal to (pi, i.e., {<Pj\(pi) = 0. 

Using a Langrangian multiplier for this constraint, the 
minimization of Ei{N ± si) with respect to ipj (while (pi 
is kept fixed) yields 

#lda| <Pj ) - Pj | <pi) = £j | <Pj) (4) 

with pj = {<pi |Hlda|</U) for the spin-up channel and 
,6j = 0 for the spin-down channel. The conventional 
conjugate-gradient formalism can be used to solve the 
above equation and obtain the minimum Ei{N ± si). 

Note, when si = 0, Eq.(3) returns to the conventional 
DFT expression. A few s/ values need to be calculated 
to obtain A/. Using Eq.©, the modified eigen energy of 
the LDA eigen state ipi can be expressed as: 

U = e^ DA + ^2h\{(PMi)? (5) 

1 

with eb DA being the original LDA eigen energy. Note, 
one could solve Eq.0 self-consistently for ipi. However, 
the self-consistent effect for bulk eigen states is rather 
small. Here we just use the original LDA wave function 
ipi, and take the expectation value of Eq.0. 

Using the method above, we have calculated 27 semi¬ 
conductor compounds (Supplemental Material), includ¬ 
ing conventional semiconductors and oxides with exper¬ 
imental band gaps ranging from 0.2 to 8 eV, covering a 
wide range of physical situations and application inter¬ 
ests. The Wannier-corrected band gaps along with the 
LDA band gaps are plotted in Fig. [I] versus the experi¬ 
mental band gaps [33(. Our LDA band gaps (Table SI) 
agree well with previous published results. They signifi¬ 
cantly underestimate the experimental values (even neg¬ 
ative for some compounds), and there is no simple cor¬ 
relation between the LDA and experimental band gaps. 

The Wannier-corrected band gaps are in good agreement 
with experiments, and the errors are on par with more 
expensive methods such as the GW Q. In the following, 
we discuss in more details for a few examples. 

GaAs is a representative of main-group semiconduc¬ 
tors, it has a direct band gap of 1.43 eV in exper- 
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iments. The LDA band gap is 0.5 eV after consid¬ 
ering the SOC. As stated before, we construct one s 
and three p projected WFs at the As site for the va¬ 
lence bands. The VBM state is purely contributed from 
the As -p projected WFs, i.e. |((/as,s|V’vbm )| 2 = 0 and 
J2 P K^As.pI^vbm)! 2 = 1, so the VBM energy correction is 
Aas , p = —0.58 eV. The CBM energy correction is purely 
from the Ga-s projected WF, i.e., |{^>Ga,s|^CBM)| 2 = 1, 
with AGa,s = 0.25 eV. Therefore, the Wannier-corrected 
band gap is 0.5 + 0.58 + 0.25 = 1.33 eV, which agrees well 
with experiments. At the bottom of the valence bands, 
the Tgu (Ti„ without SOC) state is purely contributed 
from the As-s projected WF, and the energy correction 
is Aas.s = —0.88 eV, which is 0.3 eV more than the cor¬ 
rection to the VBM. The LDA bottom valence band is 
12.84 below the VBM, so the Wannier-corrected value is 
13.14 eV, which agrees excellently with the experimental 
value of 13.1 eV 351. 

ZnO is a prototypical metal oxide with a direct band 
gap of 3.4 eV in experiments. However, its LDA band 
gap is only 0.66 eV. There are controversies for whether 
the GW can reproduce the experimental results 36, 371. 
The typical hybrid functional also gives a too small band 
gap (2.5 eV) [38}]. For the valence bands, we constructed 
five Zn -d projected, one O-s projected, and three O-p 
projected WFs. Due to the strong p-d hybridizations in 
ZnO, both the O-p and Zn-d projected WFs contribute to 
the VBM, with the energy corrections of -1.1 eV and -0.9 
eV, respectively. For the CBM, similar to the GaAs case, 
only the Zn-s projected WF contributes, and the energy 
correction is 0.75 eV. As a result, the Wannier-corrected 
band gap is 3.41 eV, in excellent agreement with experi¬ 
ments. We note that a major part of the large correction 
comes from the d component in the VBM. For II-VI sys¬ 
tems with weaker p-d hybridizations, e.g. ZnS, the d 
projected WFs only contribute ~0.2 eV to the VBM cor¬ 
rection and the LDA band gap errors tend to be smaller 
than that of ZnO. 

Another interesting quantity in ZnO is the energy po¬ 
sition of the Zn 3 d bands. The density of states (DOS) 
calculated by the LDA and the new method along with 
the experimental results are plotted in Fig. [2] The_exper- 
imental Zn 3d peak is ~7.5 eV below the VBM j39|. The 
LDA peaks are ~2 eV higher compared to experiments. 
The Wannier-corrected peak is ~7 eV below the VBM, 
which agrees with the GW results and is better than the 
hybrid functional results [38|, |4Cj. Furthermore, in ex¬ 
periments, Zn 3d is a single peak, while LDA yields two 
peaks and the low-energy peak couples strongly with the 
O 2 p states [H|. Our Wannier-corrected method reduces 
the p-d repulsion, and merges the two peaks into a single 
one in agreement with experiments. 

The above discussion indicates that the Wanner- 
corrected method not only improves band gaps, but also 
improves band energies inside valence bands. The same 
is true for conduction bands. Table Q] shows that the 


conduction band energies at the T, X, and L points are 
all corrected well by our method, for two most widely- 
studied semiconductors: Si and GaAs. 

The above results show that the new method works 
well for solids. One remaining question is whether it also 
works for molecules. In molecule calculations, we use 
the open-boundary condition for the Poisson equation 
to avoid image interactions j29| . For atoms or simple 
molecules such as LiCl, the projected WFs are just the 
eigen orbitals, so our method behaves the same as the 
ASCF method and agrees well with experiments 24|. 

Next we discuss a more complex molecule serial: poly¬ 
cyclic aromatic hydrocarbons, shown in Fig. (3[a). As 
the number of benzene-ring increases, the molecules will 
eventually become a ID system. We symmetrically pick 
half of the carbon and hydrogen atoms [circled by red 
in Fig. [3](a)] , and construct C-s, C -p, and H-s projected 
WFs on these atom sites. Other choices have been tested 
and the results are similar [to keep the symmetry, the 
WFs at the atoms uncircled in Fig. [3] (a) can also be in¬ 
cluded in Eq.([2]), with an overall factor of 1/2 applied to 
all A/]. In these cases, the WFs are no longer the eigen 
orbitals. The VBM of these molecules are purely con¬ 
tributed from the the C -p z projected WFs constructed 
by the occupied orbitals, and the CBM are purely con¬ 
tributed from the C -p z projected WFs constructed by 
the unoccupied orbitals. Figure [3] plots the experimental 
VBM and CBM [4lJ, the LDA eigen energies, the ASCF 
energy [E(N) — E(N—1) and E(N+1) — E(N)}, and our 
Wannier-corrected results. The LDA eigen energies are 
higher for the VBM and lower for the CBM by 2~3 eV 
compared to the experiments. The Wannier-corrected re¬ 
sults agree excellently with the experiments for the CBM, 
and are only ^0.4 eV lower than the experimental val¬ 
ues for the VBM. Although the ASCF energies also agree 
well with the experiments (note for benzene and naphtha¬ 
lene, the CBMs are above vacuum, so E(N + 1) cannot 
be calculated), the energy trend of the new method is 
much better. As the number of benzene-ring increases, 
the ASCF energy errors (compared to experiments) for 
both the CBM and VBM keep increasing, and eventu¬ 
ally the ASCF energies will converge to the LDA eigen 
energies for the extended ID system (the increase of the 
error is nevertheless slow, which explains why the ASCF 
works for moderate-sized molecules). However, the er¬ 
ror of our method does not increases with the molecule 
size. It indicates that the new method works from small 
molecules all the way to extended solids, and thus it is a 
general approach applicable to all systems for electronic 
structure calculations. 

In summary, we have presented a general and efficient 
method to calculate electronic structures from molecules 
to solids. This method is based on an extended Koop- 
mans’ theorem, and uses the WFs as localized orbitals 
for the fractional electron addition/removal. It does not 
have any adjustable parameters and is computationally 
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much cheaper than the hybrid functionals or GW calcu¬ 
lations. For solids, it yields accurate eigen energies not 
only for band-edge states but also for other states insides 
the bands; for molecules, it yields better trends than the 
traditional ASCF method. 
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sources of the National Energy Research Scientific Com¬ 
puting Center (NERSC) at the LBNL. Wang is also sup¬ 
ported through the Theory of Materials Project at the 
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TABLE I: The band gaps (eV) measured from the VBM to 
the conduction bands at the T, X, and L points of Si and 
GaAs. 


Compound Conduction band 

LDA Wannier-corrected Experiment 


T 15c 

2.58 

3.21 

3.34 


to 

O 

3.53 

4.09 

4.15 

Si 

X lc 

0.63 

1.2 

1.13 


Lie 

1.58 

2.17 

2.04 


L 3 C 

3.33 

3.96 

3.91 


Tic 

0.5 

1.33 

1.43 


Tisc 

3.82 

4.67 

4.72 

GaAs“ 

Xlc 

1.32 

2.14 

2.18 


X 3c 

1.52 

2.33 

2.58 


Lie 

0.90 

1.73 

1.85 


“In Ref, t.'ioi . the single group notation was used, so we average 
the SO-splitted energies of the T 15 c state. 



Experiments (eV) 


FIG. 1: The calculated LDA and Wannier-corrected band 
gaps versus the experimental band gaps [§3| for the 27 solids. 
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Energy (eV) 


FIG. 2: The DOS of ZnO from the LDA and Wannier- 
corrected calculations. The dashed line is the experimental 
UPS curve taken from Ref. j42f. Because the cross section is 
not considered in DOS calculations, the height of the experi¬ 
mental curve should not be compared with the calculations. 
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